Dropping cold quantum gases on Earth over long times and large distances 



O 
O 
(N 

> 

o 

O 

m 
(N 



X3 
O 

o 



> 
cn 
o 

o 



o 



X 
S3 



G. NandiE R. Walser, E. Kajari, and W. P. Schleich 

Institut fiir Quantenphysik, Universitat Ulm, D-89069 Ulm, Germany 
(Dated: February 6, 2008) 

We describe the non-relativistic time evolution of an ultra-cold degenerate quantum gas (bosons/fermions) 
falling in Earth's gravity during long times (10 sec) and over large distances (100 m). This models a drop 
tower experiment that is currently performed by the QUANTUS collaboration at ZARM (Bremen, Germany). 
Starting from the classical mechanics of the drop capsule and a single particle trapped within, we develop the 
quantum field theoretical description for this experimental situation in an inertial frame, the corotating frame of 
the Earth, as well as the comoving frame of the drop capsule. Suitable transformations eliminate non-inertial 
forces, provided all external potentials (trap, gravity) can be approximated with a second order Taylor expansion 
around the instantaneous trap center. This is an excellent assumption and the harmonic potential theorem applies. 
As an application, we study the quantum dynamics of a cigar-shaped Bose-Einstein condensate in the Gross- 
Pitaevskii mean-field approximation. Due to the instantaneous transformation to the rest-frame of the superfluid 
wave packet, the long-distance drop (100m) can be studied easily on a numerical grid. 
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I. INTRODUCTION 

Dropping toys on the floor is one of tiie earliest childhood 
experiences and remains a source of endless joy. Modern 
physics, as we know it today, is also based on the consequent 
pursuit of this naive amazement about the gravitational attrac- 
tion between material bodies. The contributions of Galilei, 
Newton and Einstein to the understanding of the free fall 
have fundamentally changed flie way we understand modern 
physics. 

Nowadays, the gravitational field of the Earth is under more 
intense scrutiny than ever before. One research branch fo- 
cuses on the complex geodynamics of the classical gravita- 
tional field of the Earth. From tide movements of the oceans 
and the atmospheric mass flows to minute wobbling of the 
instantaneous Earth rotation axes due to liquid core motion, 
all such effects become measurable. This is either done with 
ground based gravitometers that measure the time-dependent 
local acceleration g \ l]. Earth's rotation 12^^, 4, 5^^ 71, or in 
space where satellite-based geodesic measurements llaj^ fioll 
put tighter limits on the higher multipole moments of the grav- 
itational potential 111). Another research direction focuses 
more on the fundamental aspects of gravity which follows 
from general relativity, for example the current measurement 
of the gravito-magnetic field with orbiting gyroscopes by the 
"Gravity probe B" experiment 1 12]. This gravitational science 
on Earth and in space connects different branches of physics 
- on the atomic, macroscopic and cosmological level. It is 
vigorously pursued by American, Russian, European, as well 
as Chinese space agencies. 

The discovery of Bose-Einstein condensation fisl fl4ll and 
fermionic superfluidity flsl [T6l flTil in dilute atomic vapors 
has introduced new members to the family of superfluid con- 
densed matter systems. Outstanding features of atomic vapors 
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are that they exists at the lowest possible temperatures, i. e., on 
the nano- and pico-Kelvin scale, and their dynamic properties 
can be engineered externally. 

These unique features are now combined in an experiment 
to study Bose-Einstein condensates in /i-gravity. At the drop 
tower facility of ZARM (Center of Applied Space Technol- 
ogy and Microgravity, University of Bremen, Germany), the 
QUANTUS collaboration 1 18] focuses cuiTently on the imple- 
mentation of a fully selfcontained ^^Rb-BEC experiment that 
fits into a small drop capsule and falls repeatedly over a dis- 
tance of 100 m. This is performed inside an evacuated drop 
tower tube and results in a residual gravitational acceleration 
of Ag/g_=^ 10^^. A status report of the experiment is given 
inRef. 01. 

The physics that can be explored with the falling BEC nat- 
urally splits into two topics. First, one can consider the exten- 
sion of atom interferometry |23|2ll|23|23|] for high-precision 
inertial and rotation sensing. Due to the long unperturbed 
free-fall times of up to 10 seconds, it can be expected that 
precision of measurements can be improved accordingly. Sec- 
ond, fundamental questions regarding the quantum nature of 
degenerate gases can be studied. Due to the possibility of de- 
creasing the trapping potential significantly in a micro-gravity 
environment without the need of external levitation fields to 
compensate for gravity, lower ground state energies than ever 
should be accessible and the pico-Kelvin physics could be en- 
tered ]24]. In the resulting ultra-large condensates (10mm), it 
is possible to gain absolute control of the macroscopic matter- 
wave as optical readout and manipulation can be performed 
with very high relative spatial resolution. Furthermore, a long- 
time unconstraint expansion of a BEC allows for a measure- 
ment of the macroscopic wave function, its higher coiTela- 
tions ] ^ E6ll and probes the very concept of long-range or- 
der Exact symmetries of the system such as the Kohn mode 
iEtHSSEqH or the breathing mode ]30] can be studied in un- 
precedented detail. 

The theoretical description of any of the previously men- 
tioned topics requires a detailed modeling of the launch proce- 
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dure, the subsequent motion of the drop capsule, the dynamics 
of the comoving condensate and its time-dependent trapping 
geometry. This requires the numerical solution of the Gross- 
Pitaevskii equation for the semi-classical field amplitude. If 
we describe it on a numerical grid, which rests in the comov- 
ing frame of the BEC wave packet, then this task becomes 
most simple. Obviously, there is the need to describe the map- 
ping between this particular non-inertial frame and the other 
possible reference frames used for observation. These consid- 
erations can be extended to the quantum depletion and the pair 
correlation function IsHls^l . 

This article is organized as follows: In Section II we shortly 
review the classical physics of a single harmonically trapped 
particle falling within the drop capsule. Section III is ded- 
icated to the quantum mechanical description of a many- 
particle system, whether bosons or fermions, trapped in a har- 
monic potential and falling within a capsule in the gravita- 
tional field. We consider the unitary representations of the 
coordinate transformations on many-particle Fock space. We 
derive a harmonic potential theorem |28] and obtain the many- 
particle quantum evolution in the reference frame of the cap- 
sule, which decouples from the free-fall motion according to 
the equivalence principle. As an application, we study the 
quantum dynamics of a Bose-Einstein condensate (BEC) in 
the Gross-Pitaevskii (GP) mean-field approximation. Due to 
the instantaneous transformation to the rest-frame of the cap- 
sule, the long-distance drop (100m) can be studied easily on 
a numerical grid. In the previous sections, we have deliber- 
ately disregarded the rotation of the Earth. This is rectified in 
Section IV. In there, we obtain the main result of the article, 
i.e. the complete quantum mechanical description of the drop 
tower experiment, as well as the transformation rules for ob- 
servables between the tower- and the drop capsule frame. In 
Section V, we extend the theory for structureless particles to 
two-level atoms, which can be coupled via an external, off- 
resonant laser In the special case of equal scattering lengths 
(e.g. *^Rb), we show that the internal dynamics decouples 
from the external motion. 



II. CLASSICAL PHYSICS OF A PARTICLE FALLING 
WITHIN THE DROP CAPSULE 

In this section, we review the classical physics of a single 
harmonically trapped particle in a free-fall experiment. There- 
fore, we need to understand the classical motion of the drop 
capsule first. 

A. The drop capsule 

1. Experimental configuration 

The drop capsule is a cylindrical container that houses the 
complete setup studied in the release experiment |33|. It is 
shaped like a projectile, with a diameter of 80 cm, a height of 
2.4 m and a gross mass M. of up to 500 kg. Either, it is lifted 
110 m to the top of the depressurized tower tube (< 10 Pa) 



and released to fall freely (4.74 s), or it can be shot up from 
the ground with a powerful air-gun driven catapult, thus dou- 
bling the time available for ballistic motion. On impact, it is 
decelerated smoothly with a container full of polystyrene pel- 
lets. The same experiment can be repeated up to three times 
per day. 

Irrespective of which specific launch procedure is chosen, 
it is clearly necessary to briefly review the basic Newtonian 
physics of the drop capsule. In the simplest scenario, we may 
safely disregard any tumbling micro-motion of the capsule 
along its flight path, causing gyro-mechanical effects. Thus, 
we will model the capsule solely by its center-of-mass coordi- 
nate p = X]i=i ^i- assume that a Cartesian coordinate 
system, fixed to the center of the Earth, is a good inertial ref- 
erence frame with basis vectors denoted by {ei, 62, 63}. This 
configuration is shown in Fig. ^ For the moment, we will 
also "freeze" the rotation of the Earth and consider those ef- 
fects explicitly in Sec. II VI 

In general, the gravitational potential of the Earth Vg is not 
circular-symmetric, nor stationary. From a geophysical view- 
point. Earth resembles a drop of a viscous fluid. During the 
past eons, it evolved into an oblate ellipsoid due to its rotation. 
Consequently, the local gravitational acceleration g is orthog- 
onal to the surface of the Earth, but does not point towards the 
center, except for the equator and the poles ll34ll . Even today, 
the geodynamic activity has not subsided but remains notice- 
able in the form of tidal oscillations and liquid core wobble. 
Thus, for the purpose of modeling high-precision drop exper- 
iments, one needs to consider the general expression for the 
gravitational potential 

Vg{p,t) = -G / ' / , (1) 

Je |p-x| 

where G = 6.6742 x 10^^^ m"^ s^^kg^^ represents the grav- 
itational constant 1 35] and denotes the mass density of the 
Earth. Currently, experimental data in the form of a multipole 
expansion to the 360*'* degree is available LI 1.1 and more satel- 
lite based geodetic measurements are on the way isl 1^ fioll . 
Despite this remarkable precision in the gravitometric data, 
it remains nevertheless true that the monopole is the domi- 
nant contribution to the gravitational potential. It is propor- 
tional to the standard gravitational constant, i. e. the product 
of the gravitational constant and the total mass of the Earth 
G /^d^a;m5(x,i) 3.986 X 10^ km^s-2. 

2. Hamiltonian dynamics 

In considering the mass M of the drop capsule, it is clear 
that classical mechanics rules its dynamics. The succinct for- 
mulation of analytical mechanics follows from Hamilton's 
principle, which demands that the trajectory defined by po- 
sition p and orientation of the capsule n, extremizes the clas- 
sical action 

5,[p,n]= f At' LMt').p{t')Mt')A. (2) 

J to 
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FIG. 1 : A harmonically trapped particle with position r and mass 
m falls towards the Earth in a drop capsule with position p and mass 
M. The center-of-mass of the drop capsule p and its major axes 
coincide also with the origin of the harmonic trap and its principal 
axes Si. 



and the Hamiltonian equations of motion represent the dy- 
namics 



P = 



dn ' 



TT = — ■ 



dHc 

dp 



(7) 



In particular, applying analytical mechanics to the drop cap- 
sule in vacuum, falling in gravity, we find 



TT 



d 



P=T7' = -M — Vgip,t), 



M 



or simply Newton's equation 



dp 



d_ 
dp 



P=-—Vgip,t). 



(8) 



(9) 



B. A single classical particle trapped in the drop capsule 

For our considerations the inertial frame of reference as 
well as the comoving frame of the drop capsule are of sig- 
nificance. Therefore, we will briefly discuss them in the fol- 
lowing sections. 



subject to appropriate boundary conditions fs^. The La- 
grangian for the drop capsule with respect to an inertial 
frame is given by 

L4p,p,h,t]^M(^^-Vgip,t)^ +^hlh. (3) 

In there, we introduced the tensor of inertia 2. The rate of 
change of n corresponds to the angular velocity or spinning 
of the capsule. As mentioned in the introduction, we will 
disregard possible torques, which arise due to weak gravity 
inhomogeneities. Thus, the internal angular momentum of 
the drop capsule is conserved In = const., in other words 
n is a cyclic variable in the Lagrangian Eq. Q. To simplify 
the following discussion, we will also drop the energy offset 
nZn/2, which has no dynamical consequences. 

The global extremum of the action Sc is attained when the 
trajectory satisfies locally the Euler-Lagrange equations 

d dLr dLr 

= — -■ (4) 

dt dp dp 

The transition from the Lagrangian to the Hamiltonian formu- 
lation is based on the introduction of the canonical momenta 

= ^=Mp. (5) 
dp 

Provided the canonical momenta can be expressed uniquely 
via velocities, then the Legendre transformation of the La- 
grangian is invertible, Det{d^Lc./dpidpj) ^ 0. It defines a 
Hamiltonian function in phase space 

Hc{p,n,t) ^ Trp{Tv) - Lc, (6) 



1. Inertial frame 

The situation considered in Fig. ^ is almost analogous to 
Newton's proverbial apple dropping in the falling elevator. 
However, in addition to the gravitational acceleration, our par- 
ticle with coordinate r experiences a linear trapping force, 
possibly time-dependent. This force is derived from a har- 
monic oscillator potential 

Vt{i,t)^]^{i®$,)-vf\t), (10) 

3 

vf\t)=Y,^Kt)'^^®'^^. (11) 

1=1 

which is tied to the center-of-mass of the drop capsule p and 
is rigidly aligned along the symmetry axes of the capsule 
{ei, 62, 63}. In general, the potential can be anisotropic with 
time-dependent trap frequencies uji{t). Both information is 

incorporated in the definition of the symmetric tensor vl {t). 

As before, the equation of motion for the falling, trapped 
particle follows straight from the Lagrangian of the drop cap- 
sule, Eq. (|3}, by adding the corresponding Lagrangian of the 
trapped particle 

L[p-, P, r, i", t] = Lc[p, p, t] + Lsp[r, r, t], (12) 

L,p[r,r,t]=m(^^-Vt{r-p{t),t)~Vg{r,t)^ . (13) 

It almost goes without saying that the back-action of the par- 
ticle on the drop capsule is negligible, since the mass of the 
drop capsule Ai ^ mis much larger than the mass m of the 
atomic particle. 
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In the inertial frame, the atomic canonical momentum p is 
identical with the kinetic momentum 



dL 



sp 



(14) 



Thus, the Hamiltonian function for the single trapped particle 
is obtained immediately as 

Hsp{r,Yi,t) = pr - Lsp 

p2 (15) 

^^ + m[Vt{v-p{t),t) + Vg{v,t)] 



and the equations of motion read as 



m 



(16) 



d 



This yields Newton's equation for the particle coordinates r 

v = -^[Vt{v^p{t),t) + Vg{v,t)]. (18) 



2. Comoving frame 

The particle can deviate only a tiny distance from the cen- 
ter of the drop capsule. It is of the order of 1 cm, hence much 
less than the radius of the Earth. Consequently, it is prudent 
to make a coordinate transformation to the comoving, accel- 
erated frame and also to introduce relative velocities 



r' = r — p(t), r' = r p 
The Lagrangian of the particle now reads 
•(f' + p)2 



(19) 



L'[r',r',t]^m 



-Vt{r',t)-Vg{r' + p,t) 



and induces a canonical momentum as 
, dL' 



P = 



sp 



dr' 



m(r' + p). 



(20) 



(21) 



Now, the Hamiltonian in the accelerated frame is given by 

H',j,{v',p',t)=p'r' -L',p 

p/2 (22) 

- f- + m[Vt (r', t) + Vg{r' + p, t)] - p'p. 
2m 



and leads to the following equations of motion 

P, (23) 



r = 



sp 



dp' m 



'sp 



dr' dr 
or Newton's equation 

d 



m—[Vt{r',t) + Vg{r' + p,t)], (24) 



Due to the smallness of the possible deviation of the trapped 
particle from the trap center, we can safely expand the grav- 
itational potential at each instant in a Taylor series along the 
trajectory of the drop capsule p{t), 



Vg{r' + p,t) = [l + r' — + -ir' + ■ ■ ] t) 

= F,(p,t) + l^W(r',p,t) + Fj2)(r',p,t) 
+5V^^\r',p,t). 

(26) 

This expansion defines the gradient field and a linear potential 
contribution. 



T/W(r',p,i)=r' vW(p,t), 



(27) 



the symmetric Hessian tensor and its complete contraction 
into a quadratic potential energy 



vf\p,t) = {^^<^^^)vMtl 

V^^\r',p,t) = \{r'®r')-vf\p,t), 

A3) 



(28) 



as well as a residual correction SVg , given in Appendix lAl 
Due to its smallness, we will tacitly disregard corrections of 
this order in all of the following considerations. 



3. Pre-release dynamics 

Prior to the release at t^, the capsule is statically attached 
to the top of the tower, i.e. p{t < to) = 0, at some hight po. 
Please note that the vanishing velocity implies that we have 
deliberately "frozen" the rotation of the Earth. This is recti- 
fied in Sec. llVI Then, we can approximate Newton's equation 
Eq. (l25t as 



^-l/t(r',po,0-v^'^(po,i), 



dr 



if\po,t) {r'-r',^{t)). 



(29) 



The equilibrium position r^^ reflects the gravitational sag of 
the particle with respect to the trap center po at some instant 
and is given by 



^f)(po,i)r' W = -vW(po,i)- 



(30) 



Presumably, its time dependence is very slow. Furthermore, 
we have gathered all quadratic potential energy contributions 
into a renormalized harmonic trapping potential as 



vf\p,t) = vf\t)+vf\p,t), 



(31) 



It is interesting to note that the trapping force is weakened in 

(2) 

the direction of gravity as the curvature coefficients ' are 
negative. 
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4. Post-release dynamics 

After the capsule is released, i. e., for times t > to, it falls 
essentially according to Eq. (|9}. Due to the previously men- 
tioned sag in the direction of gravity, the equilibrated atomic 
particle will be in a position r'{to) = r'^y{to) and initiate har- 
monic oscillations about the instantaneous trap center p{t) 

^' = Vt{r\p(.t),t) = -vi'\p{t),t)r'. (32) 



5. Hamiltonian formulation and canonical transformations 

In order to obtain the transformed single-particle Hamilton 
function H'gp{r' , p', t) in terms of the old Hamilton function 
Hsp{v. p, t) of Eq. M5\ . we used the coordinate transforma- 
tion of Eq. M9\ and stepped through the Lagrangian proce- 
dure. More elegantly, one can also employ the general ap- 
proach of canonical transformations |36]. In particular, we 
will choose a time-dependent generating function 

Tn 

G(r, p", i) = [r - p(i)]p" + —[s[t) + r^{t)] (33) 

that depends explicitly on the old particle coordinate and the 
new momentum variable {r, p"}. It also depends parametri- 
cally on the coordinate and momentum of the drop capsule 
{p(t), 7r(t)}. For convenience, another, yet undetermined 
variable s{t) was introduced. This does not affect the dynam- 
ics at all, but can be used to match the energy-zero level at 
some instant. The new coordinate r" and the old momentum 
p are obtained from the generating function as 

// dG , , dG I, m , . ,^ 

r" = ^=r-p(t), p^-^p" + -.(t). (34) 

Note that besides the displacement of the coordinate r, we also 
introduced a boost in momentum space, which differs from 
the transformation described by Eq. M9\ . With respect to the 
transformed frame, the new Hamiltonian reads formally as 

i/:;(r",p",i) = i/.,(r,p,t) + ^G(r,p",i). (35) 

By inserting the coordinate transformation as well as employ- 
ing the Taylor expansion of the gravitational potential, we find 
the residual harmonic Hamiltonian 

H'Jj,{r",p'\t) = H_ +ml4(r",p(t),i), (36) 

provided the motion of the capsule is externally determined 
by the solution of Eq. (|8j and s{t, to) is an action 

sit,to)= f dt'LMt'),Pit'),t']-p7r\l. (37) 

Jto 

By construction it is clear that we must recover the identical 
Eq. j32l l from Hamilton's equations in the new coordinates 



We would like to point out that these considerations are not 
purely academic but will enlighten the analogy between the 
classical treatment and the quantum mechanical description 
of the falling, trapped ultra-cold quantum gas. 

6. Order of magnitude estimates 

In order to estimate the magnitude of the Taylor coeffi- 
cients of the gravitational potential, one can evaluate them 
at the surface of the Earth p = — 6372 km and ob- 
tains = 9.81 ms-2, uf' = -1.54 x 10"*^ s^^ and 
Vg"^^ — 2.42 X 10^^^ m^^ s^-^. In there, we made the ap- 
proximation of an isotropic gravitational potential. It is worth- 
while pointing out that Vg^^ is not yet identical to the gravita- 
tional acceleration g, which also includes centrifugal correc- 
tions (see Sec. lIV> . 

The important consequence of Eq. (j^^ is that in first or- 
der Vg^\ the internal harmonic oscillator motion is decoupled 
from the center-of-mass motion of the drop capsule. Only 
when considering the weak quadratic correction Vg , one can 
observe a coupling of these motions as a result of the decrease 
of the harmonic oscillator frequency along the gradient of the 
gravitation. The following quantum mechanical calculations 
will obviously lead to the same conclusions due to Ehrenfest's 
theorem, where in case of quadratic Hamiltonians expectation 
values coincide with the classical trajectories. Only by con- 
sidering the third order correction 6Vg , we will find an addi- 
tional dynamical mixing between the classical trajectory and 
the quantum mechanical expectation value. 

III. QUANTUM PHYSICS OF IDENTICAL ATOMS 
FALLING WITHIN THE DROP CAPSULE 

So far, we have considered a single classical particle that 
is harmonically trapped and falling in the gravitational field. 
Actually, we are interested in the behavior of a freely falling, 
trapped atomic Bose-Einstein condensate (BEC), which re- 
quires a quantum field theoretical description. However, we 
do not want to constrain our discussion solely to condensed 
bosonic gases, but would rather like to treat the general situa- 
tion of an ultra-cold degenerate quantum gas, whether bosonic 
or fermionic. 



A. Quantized atomic fields 

Before we proceed to the quantum field theoretical descrip- 
tion, which is intrinsically tied to the picture of second quan- 
tization, we briefly revisit the most prominent relations of first 
quantization. In general, the commutator of position and mo- 
mentum operator f and p is 

[r,p] = ihl. (39) 



In the remainder of this article, we will omit the hats for op- 
erators in first quantization, e.g. write r and p. The derived 
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commutation relation for the components of the angular mo- 
mentum 1 = r X p satisfies the standard angular momentum 
algebra 



Ik, h] = ifl^klm Ir. 



(40) 



defined through the completely anti-symmetric Levi-Civita 
tensor ekim as structure constants. Most of the time, we will 
use the position representation in this article, hence will inter- 
changeably refer to r and p = also as the operators of 
position and momentum, when acting on Hilbert-space func- 
tions. 

In order to incorporate the correct quantum statistics ac- 
cording to the Pauli principle, we use the language of second 
quantization and introduce the basic quantities in the follow- 
ing. The field operator in the position representation will be 
denoted by a(r) and it annihilates a particle at the position 
r, while its hermitian conjugate at(r) creates it there. In the 
present discussion, all operators are given in the Schrodinger 
picture unless indicated otherwise. Thus, the quantum field 
satisfies a commutation relation for bosons (— ) and an anti- 
commutation rule for fermions (+) 



in terms of single-particle Hamiltonians Hsp{r,p,t), pair- 
wise interactions Vp{^) and, with decreasing relevance, higher 
order effects, which can be neglected here The re- 

cently discovered Efimov-resonances represent a very un- 
usual exception to this rule, where genuine three-body effects 
become important. However, the principle of translational in- 
variance can still be assumed for those potentials. For the 
present discussion, we will disregard such peculiarities, how- 
ever Obviously, we want to assume in here that the quan- 
tum dynamics of the A^-particle state follows the A^-particle 
Schrodinger equation 



(47) 



Now, if this is translated to the view point of second quan- 
tization, one obtains the standard Hamilton operator 



Hsp{t)^ j dVat(r)i/,p(r,p,i)a(r), 
Vp^\ j d\a\vi)a){T2)Vp{vi~V2)a{v2)a{vi). 



(48) 



[a(r),at(r')] -^(r-r'). 



(41) 



For commutators, we will drop the " in the following. 

Now, we can proceed to single-particle operators in Fock 
space representing the center-of-mass coordinate R, the linear 
momentum P, as well as the angular momentum L and the 
total atom number N as 

R = y" d^rat(r)ra(r), P = ^ dVa^(r)pa(r) 



In low-energy field theory, the pair potential Vp is often mod- 
eled not by a finite range potential, but it is assumed to be 
proportional to a delta function. This is called the contact- 
or pseudo-potential approximation. It is very useful in cer- 
tain situations, but needs careful consideration when going to 
higher order calculations Issil . For the present discussion, we 
can postpone this question and simply demand translational 
invariance. 

In principle, one can represent any state in Fock space as 



f aj f ^ u ^ ^^^^ |*) = / d'^^Hri, • ■ • , tat) (n, . . . , | vI/(^)), 

L / d-'raT(r)la(r), (Tr a' (r) a{r). 



(49) 



Using the commutation or anti-commutation relation of 
Eq. (I41> . one can easily verify that those operators form a 
closed algebra 



|ri,. 



R,P 



R,A^ 



P,iV 



R, r2L 



i?i r2 X R, 



ihN, 
P.fJL 



(43) 

= 0, (44) 
ifiVlxP. (45) 



It is indicated by Eq. (I44> that this algebra of Euclidean trans- 
formations can be represented in the A^-particle sector and 
O represents an arbitrary rotational axis in three-dimensional 
Euclidean space. 

In typical BEC experiments, as well as in recent experi- 
ments involving degenerate fermions, the particle densities are 
of the order 10^^ cm^'^. Thus, we can express the energy of a 
dilute atomic gas of N particles in form of a cluster expansion 



JV 



N 



(46) 



as a superposition of different A^-particle states |ri, . . . , tn) 
created out of the vacuum \v) and weighted by the prop- 
erly symmetrized or anti-symmetrized -particle amplitudes 
vI/W(ri,...,r^) = (ri,...,rAr|*W). 

Proceeding formally, it is also well-known that the 
complex-valued Lagrangian for Schrodinger fields L3^ l40l 

EH 

C[^it),i'it),t] = {^it)\zh--H{t)\M> it) ) (50) 
leads via a Euler-Lagrange variational procedure 

^{t) = ^^^h^itr, m^^, (51) 

also to the Schrodinger equation in Fock space 

^n^^\^{t))=Hit) (52) 

finally. With this fundamental equation of motion, we have 
now established all the basic concepts and can turn to the 
transformation of a wave function from an inertial frame of 
reference to an accelerated, comoving frame of reference. 
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B. Static representation of states in Euclidean frames of 
reference 

In this subsection, we will briefly review basic transforma- 
tion theory for wave functions under spatial translations and 
rotations of the coordinate system ll4H43ll44ll45ll46il . 



1. Three-dimensional Euclidean space 

The basic operation of a translation in space is given by 

Tn[A=-r + Tl- (53) 

The whole set of translational operators forms an Abelian 
group i.e., T-n^T-ji^ = T-ji^^-ji.^- Each element of this con- 
tinuous Lie group is indexed by a three-dimensional vector 

n. 

A rotation around an axis for a time t is characterized 
by a total rotational vector Q = i. This operation will be 
denoted by Mq. The set of all those isometric, linear trans- 
formations forms the infinite-dimensional Lie group SO (3), 
which excludes the possibility of reflections. The entirety of 
operations establishes the Euclidean group and a particular el- 
ement is denoted by 



r' = TuWqv] = MQV + n. 



(54) 



If we represent the basic position vector in Cartesian coordi- 
nates as r = f^^i7 then one finds a three-dimensional 
orthogonal rotational matrix 



Mqy 



Q X r 



(55) 



in terms of the generators of the rotation {^1,^2,^3}- For 
convenience, we have introduced those 3 x 3-matrices as 
complex- valued angular momenta in units of Ti, which satisfy 
the angular momentum algebra of Eq. ( I40l l. No quantum as- 
pects are implied by this choice of notation (see also Appendix 

m. 



2. Single-particle Hilbert space 

The action of the translational operator in single-particle 
Hilbert space 

Un = e-'/'^'^P, (56) 
can be seen most clearly by applying it to position eigenstates 
|r)' = C/K |r) = |TK[r])==|r + 7^). (57) 
Obviously, the operation 

Uv = e:'^'^\ (58) 

is just a boost in momentum space, when acting on momen- 
tum eigenstates, i. e.. 



The unitary representation of the rotation 

Uq = e-'l^ ^\ (60) 

is induced by the angular momentum operator 1 of Eq. (|^). 
When acting on position as well as momentum eigenstates, 
one gets 

|r)' = [/Q |r) = |MQr), (61) 

|P>'-C/S Ip) = |MqP)- (62) 

Thus, a representation of a general element of the Euclidean 
group in Hilbert space is 

U^UsUv Un Uq, Us^e'/"^. (63) 

This definition of a general group element is not unique, but 
admits the inclusion of an arbitrary phase factor Us, such that 
the transformed position state reads 

|r)' = U\r)^ e'l'^ [5+P(A/cr+i^)] j^j^^ ^ (54) 

Such a ray-representation will become physically important 
when considering the dynamic evolution of a quantum state, 
and the accumulated phase will reflect the classical action. 



3. Many-particle Fock space 

We will now proceed to a wave-function \^) in many- 
particle Fock space and transform it to the wave-function 
in another Euclidean coordinate system 

m^Um'. (65) 
A general element of the Euclidean group in Fock space 



is now given by the following representation 



(66) 



(67) 



|p)' = C/p |p) = |Tp[p]) = Ip + T') 



(59) 



where the single -particle operators R, P, L and N have been 
introduced in Eq. ( I42> . 

In the following context, we will not need the action of op- 
erators U on state vectors I^P) in Fock space, but rather on the 
scalar amplitudes of the quantum field in the position repre- 
sentation a(r). This is straight forward using the commutator 
relations of AppendixlHlin case of the phase operator Ug and 
momentum boost operator JJ^ 

Ula{v)lJg^Usa{v), ul d{r) U^:, ^ U-p air). (68) 

For the spatial shift operator Uj^ as well as the rotational op- 
erator Uq, we need to recall that we are dealing with a scalar 
field. Then one finds that the unitary operations in single- 
particle Hilbert space U-ji and Uq correspond to the inverse 
operations in Euclidean space and Mq^, i. e., 

?74 a(r) C/^ = Utz a(r) = aiT^^h) = a(r - 71), (69) 
f/^ a (r ) ;7q - C/g a (r) = a ( Afg 1 r ) . (70) 
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C. Dynamics in comoving frames of reference 

I. Frame transformation from the inertial to the comoving frame 

The Euclidean coordinate transformations of the previ- 
ous section where determined by ten static parameters, 
{S, TZ, V, Q}. By considering an arbitrary time-dependence 
of these parameters, we can introduce a canonical frame trans- 
formation that includes the Galilean transformation as a spe- 
cial case 142). Currently, we will focus on non-rotating frames 
iQ — 0) for simplicity, i. e.. 



U{t) — Us(t) U-pf^t) Un{t) 



(71) 



and lift this restriction in Sec. llVI 

Our goal is to determine an equation of motion for this co- 
moving, accelerated frame such that the residual motion of 
the atomic gas within the frame is free of any non-inertial 
forces. Given the Schrodinger equation (I52> holds in the in- 
ertial Earth-centered frame, then we find another realization 
l^*) — U (t) in an accelerated frame of reference 



in-\^{t)y^H'{t)\^{t)y, 



(72) 



with the transformed Hamiltonian 



H\t) =. U\t) i-ih- + H,p{t) + Vp) U{t). (73) 

In here, the first contribution to the Hamiltonian is a gauge po- 
tential, in analogy to the fictitious forces that would appear in 
the classical Hamiltonian function in terms of the time deriva- 
tive of the generating function, cf. Eq. ( I35> . It can be evalu- 
ated using Eqs. jB13llBT^ and it generates additional gauge 
forces 

uHt){-ih^)u it) = is + vn)N + rfi-iiP. (74) 

at 

The second contribution to the Hamiltonian implements the 
spatial translation and momentum boost 

U\t)Hsp{t)U{t) (75) 
= j dV a){v')H,p{v' + n,p' + V,t)a{v'), 

where it is implicitly assumed that the one-particle Hamilto- 
nian Hsp{r, p, t) can be expanded into a power series in r and 
p, respectively. 

The third term of the Hamiltonian Eq. ( I73> is not affected 
by the frame transformation. 



(76) 



as the interaction of Eq. J48t is hermitian (number conserv- 
ing), a local operator, and translational invariant by construc- 
tion in Fock-, as well as in two-particle Hilbert space. 



This property is also the basis of Kohn's theorem frf], i. e., the 
exact decoupling of the center-of-mass motion and the relative 
excitations in maximally quadratic single-particle Hamiltoni- 
ans. It is also known as the harmonic potential theorem |2§] 
and has been used in the context of cold dilute quantum gases 

mill. 

Combining these results, we find the transformed Hamilto- 
nian in the accelerated frame 

H'{t)= J dVat(r')i7^p(r',p',i)a(r') + 14, (78) 

and all gauge contributions are contained in the single-particle 
Hamiltonian 



H',p{r',p',t) ^Hsp{r' + n,p' + r,t) 

+ s + vn + rr' -Tip'. 



(79) 



2. Application of the frame transformation in the harmonic 
approximation 

We will now apply the preceeding considerations to the 
falling trapped interacting many-particle system in the drop 
tower The single-particle Hamilton operator follows from the 
classical considerations of Eq. il5i and is given by 

Hsp{r, p, t) = ^ + m [Vt (r - pit), t) + Vg{r, t)]. (80) 

By splitting the transformed Hamiltonian Eq. il9\ into con- 
stant, linear and higher order polynomial contributions in 
terms of the position r' and the momentum p', one finds 

H'ir', p', t) = (ap' + /3r' + 7) + ^ + mVt (r', t) 



2m 



Vg{r' + TZ,t)- Vg{n, t) - V}'\r\ 7^, t) 



(81) 



with coefficients 
V 



n, (3 = v 



m—[Vt{n~p,t) + Vg{n,t)], 
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2m 



m [Vt{n-p,t) + Vg{n,t)] + s + rn. 



(82) 

If the constraints are satisfied identically, i. e., {a, f3, 7} 0, 

one obtains the single-particle Hamiltonian, if we disregard 

(3) 

again terms of third and higher orders SVg in the accelerated 
frame. 



Hip{r',p\t)=^ 
^ 2m 



•,Vt{r\nit),t) 



(83) 



provided the frame parameters as well as the center-of-mass 
of the drop capsule {TZ{t) , V (t) , p{t)} satisfy the classical 
equations of motion 



[N,Vp] = [R,Vp] = [P,Vp] = 0. 



(77) 



jy _ dLsp 

dil 



Mp, 



mlZ, V 



dp' 

dLgp 

on ■ 



(84) 
(85) 
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We can express the latter line as Newton's equation of motion 
for TZ, 



d 



[Vt{n-p,t) + Vg{n,t)]. 



(86) 



With Eqs. (I84> and (I85> . we have recovered the Lagrangian 
L[p, p, TZ, TZ, t] of Eq. (I12t . and one can express the accrued 
dynamical phase as an action integral 



sit,to)= / dt'Lsp[n,n,t']-v7i\ 



(87) 



like in Eq. (|2j, where we started from, cf. also Eq. i37\ . 

The harmonic trapping potential Vt that appears only in 
the residual Hamiltonian, Eq. ( I83> . has been modified by 
quadratic corrections of gravity and was defined in Eq. ( I31l l. 
The third order correction to the gravitational potential SVg^"^ 
is minuscule. On the time scale of the free fall experiment, 
one may safely disregard it. Thus, we are left with a quadratic, 
time-dependent Hamilton operator and the harmonic potential 
theorem applies |281. 

As before, the motion of the center-of-mass coordinate TZ, 
Eq. il, can be understood much more clearly, if we expand 
the gravitational potential around the center-of-mass of the 
drop capsule p. For the deviation 6TZ{t) — 'TZ{t) — p{t) after 
the release of the capsule, one finds a pure harmonic oscilla- 
tion 



(57t(i) 



{p{t),t)5n(t). 



(88) 



D. Ehrenfest's theorem 



At this point, we would like to confirm the physical inter- 
pretation of the frame parameters ??.(<) and 7-*(<) as center-of- 
mass and momentum coordinates, respectively. Clearly, they 
have been introduced as such with dimensions of length and 
momentum. Using Ehrenfest's theorem, we are able to verify 
this. 

The temporal evolution of the expectation value of any 
time-independent operator^, i. e., {A) = {^{t) \ A\^{t)) fol- 
lows from the Schrodinger Eq. ( I52> and gives 



ih-{A) = {[A,H(t)])- 



(89) 



By using the bosonic or fermionic commutation relation, 
Eq. (14 1> . we can work out the commutators in the Schrodinger 
picture 



p 

m ' 



l[H{t),V]^-v^^\p{t))N 



(90) 
(91) 



,-,(2) 



{p{t),t) (R-p(i)iV), 



where we have discarded third order corrections. In analogy 
to Eq. ( I88> . we can define a deviation from the drop capsule 



coordinate as ^(R) 
according to 



(R) - p{t)N and find that it evolves 



= -f,(2) 



ipit),t)d{R) 



(P)=m(R), (92) 



where we have introduced the particle number N — (N). In- 
deed, Eqs. ( I32t . ( I88> and (|92} are identical in form, which 
is just Ehrenfest's theorem. Furthermore, if we assign the 
identical initial conditions to the frame parameters TZ{tQ) = 
(R{to))/N and V{to) = {P{to))/N, then we can identify 
them with the center-of-mass coordinate and total momentum 
of the atomic ensemble. 



E. Application of the comoving frame transformation to 
mean-field theory 

All the calculations that were presented so far are generally 
valid on the many-particle level using the language of sec- 
ond quantization. We have seen that the underlying physics of 
frame transformations in translational invariant systems rests 
exclusively on the single-particle level. If we furthermore re- 
strict our scope to harmonic potentials, everything can be un- 
derstood in terms of classical physics, in principle. 



1. From quantum fields to classical fields 

For the practical purpose of studying the internal dynamics 
of Bose-Einstein condensates or superfluid fermionic gases 
in //-gravity, it is usually sufficient to restrict oneself to the 
mean-field picture. However, the application of the presented 
frame transformations to an extended mean-field approach in 
the case of bosons |31] or fermions |48] is straight forward. 

In the case of a bosonic gas, the Gross-Pitaevskii (GP) 
equation provides a very successful description of many dy- 
namical and static properties of the BEC |49]. To derive the 
GP equation in the accelerated frame, we start from Heisen- 
berg's equation for an operator A, which is time-independent 
in the Schrodinger picture 



in-AH{r,t) = [AH{r,t),H'Hit)], 



(93) 



which can be obtained from the Schrodinger equation in 
the accelerated frame Eq. (^2} and the Hamilton operator 
Eq. ( I78> . The representation of operators in Heisenberg pic- 
ture requires the introduction of the time-ordered (T) evolu- 
tion operator 



ijf(r,t)=ft(t,0)i(r)f(i,0), 
-i/h / dTff'(r) 

"'0 



T{t,0) = Texp 



(94) 
(95) 



In particular, one finds for the field operator 



aHiv' ,t). 



(96) 
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In the case of a dilute weakly correlated BEC, it is usually 
sufficient to stay within the mean-field approximation. This 
is the classical limit and replaces the quantum field a//(^, t) 
by a complex valued field a{^,t). This approximation is also 
consistent with the use of a pseudo potential for the pair inter- 
action i. e., 

— -m = ^^m- (97) 

m 

It uses the s-wave scattering length in vacuo and basically 
summarizes all the binary scattering contributions to the Born 
series that are still present in the full Heisenberg Eq. ( I96> . 
With these assumptions, one arrives at the GP equation in the 
accelerated frame 

indta{v',t) = [iJ^p(r',p',i) + K|a(r',t)|2] a(r',i), 

(98) 

with the Hamiltonian of Eq. ( I83> and the classical frame equa- 
tions Eq. dSSt . 

2. Numerical study of the long time evolution of a freely falling 
condensate 

In order to illustrate the benefit of the transformation to 
a comoving frame, we present a simple numerical calcula- 
tion of a BEC in a time-dependent trap after this transfor- 
mation. Within the GP mean-field picture, we calculated 
the corresponding time-evolution for two time-dependent trap 
configurations. For simplicity, we considered a quasi one- 
dimensional BEC of *^Rb atoms and solved the corresponding 
GP equation. This is possible for a very prolate trap configu- 
ration, ^{t) — LOr/i^z (t) ^ 1, with the radial and longitudinal 
oscillator angular frequency ujr and ujz{t), respectively. The 
radial motion is effectively frozen out. Then, we can separate 
the ground state of the harmonic oscillator in radial direction 
and integrate out this degree of freedom. This yields 

,ndta = (-—dl + a^^^^z^ + ^'\a\A a. (99) 
\ 2m 2 / 

In there, a = a{z, t) is the order parameter within the one- 
dimensional GP equation in the comoving frame of reference. 
A constant energy-offset arising from the integration over the 
radial part has been removed by switching to an interaction 
picture, k' is the coupling constant modified due to the ra- 
dial integration. In Figs. Hand |3l we show the time evolu- 
tion of a *^Rb condensate consisting of = 10000 atoms 
for two situations: In the first case, the free time evolution 
is studied, i.e. the trap is switched off instantaneously at 
t = 0, when starting the drop experiment, which corresponds 
to u!z{t) — LOz Q{—t) with the well-known Heaviside func- 
tion. Due to the duration of about 5 seconds, the BEC can 
expand largely. Our chosen set of parameters is the atomic 
mass of *^Rb m = 86.9 amu (unified atomic mass units), 
u!z — 1/s, and uJr — 2tt 100/s. In the figures, we used scaled 
quantities, i.e. time is normalized to the longitudinal oscillator 
frequency, t = t uj^j^-K, lengths are measured in harmonic os- 
cillator units a^o — yjfi/muoz — 27fim. The scaled coupling 



constant is given by R' = 2j{0) ag/aho, where the s-wave 
scattering length ~ 5.8 nm. For the ^^Rb parameters see 
for example |50]. In the second case, we simulated the time 
evolution of a condensate, which is displaced from its equi- 
librium in the trap by one scaled length unit at the instant of 
release. The oscillator motion decouples from the free fall tra- 
jectory. In here, we assumed a constant trap angular frequency 

UJz = 1/S. 



500- 
400 

'"-^ 300 ■ 
200- 
100- 



-200 -100 100 200 

X 

FIG. 2: Free expansion of an initially trapped BEC over 5 seconds. 
Length is measured in harmonic oscillator units. The dashed line 
shows the density of the BEC at the instant of release, the solid line 
shows the condensate at the end of the free expansion. 




-15 -10 -5 5 10 15 
X 



FIG. 3: Oscillatory motion of a condensate that is initially displaced 
by a;o = 1 a/io from the equilibrium position. The Kohn mode os- 
cillates with the trap angular frequency ljz and decouples from the 
free fall dynamics. The solid line shows the density of the BEC at 
equilibrium, the dashed lines correspond to the turning points. 
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IV. MANY IDENTICAL ATOMS WITHIN THE DROP 
CAPSULE IN A ROTATING FRAME 

In the previous sections, we have deliberately omitted the 
rotation of the Earth from our discussion in order to simplify 
the algebra and focus on the essential physics. In particular, 
we only considered the transformation from the Earth-fixed 
inertial frame {e^}, located at its origin, to the accelerated 
frame of the freely falling, center-of-mass coordinate TZ{t) of 
the atomic cloud (see Fig.QJ. 



A. Classical physics of the drop capsule in a rotating frame of 
the Earth 

1. Classification of the three important frames of reference 

In fact, in our considerations three different frames of refer- 
ence are of significant importance. Besides the inertial frame, 
a real experiment obviously requires a description in a coor- 
dinate system that is aligned with the corotating drop tower 
tube. Moreover, in the frame comoving with the drop capsule, 
the experiment can be described in a most simple way. In 
this frame of reference, non-inertial forces can be eliminated. 
It is of special importance to describe the mapping between 
the drop tower frame and the capsule frame, since both sys- 
tems can be used for observation, in principle. In the inertial 
frame, a point-like particle is characterized by (r, p), and the 
drop capsule is at position p. In the corotating drop tower 
frame, with origin at the bottom of the tower, we indicate 
these quantities with (f p') and p' , while in the drop capsule 
frame the particle's coordinate and momentum are (r",p"). 
A schematic sketch of the three systems is depicted in FigE] 



2. Characterization of the drop tower frame 

In the rest frame of the drop tower, the coordinate system 
is aligned with the rotating tower principal axes denoted by 
{ej{t)}. This is depicted in Fig.|5] 8 denotes the angle be- 
tween the rotational axis of the Earth and the plummet at the 
geographic location (0, O) of the drop tower 5J. Due to the 
ellipsoidal shape of the Earth, the plummet does not point to 
the center of the Earth. The oblateness is modeled by the half 
axes a and h. The instantaneous principal tower axes, as well 
as the direction towards the base of the tower 5R(t) 



ei{t) = M^{t) M0e3 Mee2 e^, 
U{t)^ M^{t)M^,,Me,, 63, 



e, = n,-K e„ (100) 
U = Q.,y,lR., (101) 



are obtained from the inertial axes of the Earth {e^} by align- 
ing the §3 -axis of the tower along the direction of the plum- 
met, followed by a rotation (p to the longitude of the tower at 
ZARM/Bremen, as well as the diurnal rotation A/j (t) around 
the Earth axis Jlj = fij 63 with an angular frequency fJ, = 
7.2 X 10~^ 1/s. No further time-dependence caused by geo- 
physical effects, like tidal forces, precession, or wobbling due 
to the liquid core motion will be considered. 




rotating drop tower frame (r'.p') 



capsule frame (r",p" 




inertial frame (r,p 



FIG. 4: Schematic view of the three different reference frames as 
seen from above the North pole. In the inertial frame, located at the 
center of the Earth, a particle is characterized by (r, p), and the drop 
capsule is at position p. In the corotating drop tower frame, with 
the origin at the bottom of the tower, these quantities are denoted 
by (f',p') and p' , while in the drop capsule frame the particle's 
coordinate and momentum are (r", p"). 



Now we are in the position to represent all vectors either in 
the inertial Cartesian basis of the Earth or the corotating frame 
i. e., 

3 3 
i=i 3=1 

3 3 

by introducing the matrix representation of the orthogonal ro- 
tation M{tY M{t) = 1 defined by Eq. ( fTOOl . It is obvious 
that the components of the Earth rotation as well as the pointer 
to the base of the tower 

= (-sin 61,0, cos 61), (104) 
= SR (sin (9-61), 0, cos (9- 61)), (105) 

must be time independent in this basis. 

3. Dynamics of the drop capsule 

All distances will be measured relative to the base of the 
drop tower SR(i), i. e.. 



p = p' + U = M{t) p = M{t) (p' + SR). 



(106) 



By differentiation, one finds the relation between the veloci- 
ties with respect to the two frames as 



p = M{t) ip' + n,xp) 



(107) 
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FIG. 5: A harmonically trapped particle with position r' and mass 
m in the drop capsule with position p' and mass A4. The position in 
the capsule frame is denoted by r". The drop tower axes {ei} per- 
form a diurnal rotation around the Earth axes flf,. 9 denotes the angle 
between the rotational axis of the Earth and the plummet at the geo- 
graphic location (0, &) of the drop tower The oblate ellipsoidal 
shape of the Earth is modeled by the half axes a and b. 



In the rotating frame, the Lagrangian of Eq. (|3} for the relative 
position p' as well as the relative velocity p' reads 



L'Ap',-p',t] ^ M ( ^'p' + ^^^'p^' - v,mt)-p,t) 



= ^[^- VciP) ~ Vg{Mit)p,t)^ + n,{p X f>'M). 

(108) 

Here, we have introduced the centrifugal potential as 



(109) 



The canonical momentum tt', that is conjugate to the coor- 
dinate p', is given by 

BT ' 

ir' = -^^Mif)' + n,xp) (110) 
op' 

and implies a Hamiltonian in the new coordinates 

K{p',iv')^n'-p'~L', 

7f'2 - (111) 

= ^+MVg{Mit)p,t)-nX, 
where the angular momentum \ = p x tt' has been intro- 



duced. Via Hamilton's equations 



diT'_ 
dH' 



M 



d 



X P, 



(112) 



(113) 



we are lead directly to the well-known Newton equations in 
the rotating frame | 



p = -2n,x p' 



d 



— [V,{p) + VgiM{t)M], (114) 



where the first term is the Coriolis force, followed by the cen- 
trifugal force of Eq. ( I109> and the gravitational force as de- 
fined before in Eq. (|9j. By expanding the gravitational poten- 
tial to second order around the base of the drop tower (see 
Fig.|5jl, we obtain 



-2n, X -p' - 



(115) 



Here, we have introduced the effective gravitational accelera- 
tion and rotated Taylor coefficients 



r,(2) 



g(5i) = (SR) + fl,x {n, xU), (116) 
Af(t)TvW(5?), (117) 
(SR) = M{t)'^v^^\n)M{t), (118) 



at the surface of the Earth, g is normal to the surface and 
takes into account the ellipsoidal shape of the Earth. In order 
to obtain a simple solution of Eq. (^^Jl, one can neglect the 
centrifugal term, since it is of the order fij^ « 10~^m/s^, 
see Appendix C. An estimate of the Coriolis deviation in §2- 
direction after a drop distance of 100 m yields approximately 
2 cm, which is not negligible. 



B. A single classical particle in the drop capsule 

1. Transformation from inertial- to tower coordinates 

Now, let us get back to the harmonically trapped particle 
at r inside the drop capsule p and consider its motion with 
respect to the rotating base of the tower at analogous to 
Eq. ( II9I 1. In this case, coordinates and velocities are given by 

r = A/(t)f = M(i)(r' + SR), (119) 
r M{t) (f' + n,x r). (120) 

Note, that this description refers to the tower coordinates. 

The trap rests in the drop capsule, which is moving on the 
classical trajectory given by Eq. (^^i. From the Lagrangian 
in the inertial frame, Eq. ( I13I I. we obtain the Lagrangian for 
the particle coordinates relative to the base of the drop tower 

L;p[f',F',i] =m[i(f' + f2, xf)2 

- VtiMit){r' - p'),t) - Vg{Mit)iv' + U),t) 

(121) 
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The total Lagrangian, which includes the Lagrangians both of 
the drop capsule as well as the single trapped particle, reads 



L' = L' + L' 



sp 



(122) 



with L'^ as defined in Eq. ( I108> . The linear canonical momen- 
tum is 



di' 



m (r' + r2j X r) 



(123) 



and thus one finds for the Hamiltonian 

Hip{r',p',t) + V,mW - P'),t) 

+ Vg{M{t){f' + K), t) - p' {n, X f) . (124) 
Then Hamilton's equations of motion read 
. . dH' 



dp' 



fJ, X r, 

m 



(125) 



^' = = " 9r + + P' ^^^^^ 

correspondingly. Finally, Newton's equation is given by 

f' = - ^ {Vt{M{t){r' - p'),t) + V,{M{t){v' + 5R),t)) 
- X {n, X r) - 2n, X f'. 

(127) 

If we expand the gravitational potential up to second order, we 
obtain 

d 



f' = - 2 X F' ~ g(p) - —Vt{v' - p', p, t) 



2n,xf'-5(2)(p,i) (f'-f' (i)). 



(128) 



As in the non-rotating situation of Sec. lllB 3l we get an equi- 
librium position v'^q{t), which is defined by 

v^^Hp,t)r',qit)^~g{p)+v^^Hp,t)p' (129) 

and a modified, appropriately rotated harmonic trapping po- 
tential 

v(^^{p,t) = vi'\t)+v^^Hp,t)+vi'\ (130) 

vi^\t) = Mit)^vi^\t)Mit), (131) 

Vtitp,t) = l{^^$)-f}^^\p,t). (132) 

The equilibrium position r'^q{t) is governed by the trajectory 
of the drop capsule p', however there is a gravitational sag. 



2. Hamiltonian formulation of the canonical transformation from 
tower- to capsule coordinates 

In order to obtain the classical Hamiltonian of a single 
trapped particle in capsule coordinates, we introduce the 
canonical transformation 

G(r', p", t) = [M (f - p'm p" + [s'it) + r'n'it)] , 

(133) 



which yields 



(134) 



M 



In this particular gauge, the one-particle Hamiltonian is 
given by 



i/:;(r",p",t) ^ i?^,(f',p',t) + |G(f',p",t) 



(136) 



2 m 

provided the equations for the capsule, Eqs. and ( II 13t 

are satisfied and s'(t, Iq) is an action. 



s'(t,io)= / dt'L',[p'rp',t']-p'TT'\l. (137) 

■/to 

Obviously, this Hamiltonian is identical to Eq. ( I36t . 



C. Many identical atoms in tlie drop-capsule 

1. Frame transformations between tower- and capsule coordinates 

The calculations of Sec.|ni]excluded the effects of rotation. 
In here, we will rectify this and obtain a general frame trans- 
formation that eliminates external forces and torques from 
Schrodinger's equation. In particular, this is applied to the 
problem of gravitational acceleration and Earth's rotation. 

We have already introduced the generic many-particle 
Hamiltonian with translationally invariant, binary interactions 
in Eq. ( I48l l. However, we have not exploited the fact that only 
the s-partial wave contributes to the two-particle scattering at 
low energies. Thus, we want to assume that the relative angu- 
lar momentum is a good quantum number. This can be mod- 
eled by an interatomic potential that is only a function of the 
relative distance, V^(|^|). In the case of fermionic particles, 
this means there is no s-wave interaction due to the Pauli ex- 
clusion principle and the gas becomes ideal at low energies. 

Let us start with the full many-particle Hamiltonian 



(138) 



V ^- 
f 2 



d6r'at(rl)at(r',)K'(|rl - v'^\)a{v',)a{v',) 



in the drop tower frame. 

Now, we have to transform the Schrodinger state from the 
drop tower frame, to the rest frame of the drop capsule, 

\^{t)) = U {t) {t))' . This requires to augment the unitary 
frame transformation from Eq. ( I71> of Sec. 1111 Cl 



U{t) - Ugif^t) Uv'{t) Uii'it) Ull^ty 



(139) 



with 



;7M(t) = e-''^^'^'^^ e'~'''^">"'''^ e-'/^'^'"'^. (140) 
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In there, we performed the same rotations as in Eq. dlOOt . 
and in particular, the diurnal rotation around the Earth axis, 
Q{t) — f2j t, is accounted for. Moreover, we would like to 
point out again that we transform from the drop tower frame, 
which is corotating with the Earth, to the instantaneous rest 
frame of the drop capsule, which is not rotating. Therefore 
we need to apply Ulj(^^-^ rather than UM^t)- 

The transformation rule for the Hamilton operator was 
given in Eq. ( I73> and consists of three contributions. The first 
one contains merely the gauge contributions that arise from 
the time-dependent frame parameters 

U\t){~ih^)U{t) =0' + P' iC) N + P'iM^R.) 
at 

(141) 

For obtaining this result, we have used the basic Eqs. ( I44t . 
( IB13t and (IB14> . The second contribution is the transformed 
single-particle Hamiltonian 

UHt)H^t)U{t) = 

dV" a^r")H',j,{M^r" +n\M^p" + V\t)a{r"). 

(142) 

Finally, the third contribution is simply the transformed binary 
interaction potential. 

U^t)v;U{t) = v; = Vp. (143) 

It is left unchanged due to particle conservation, the local 
character of the potential, as well as the translational and ro- 
tational invariance, i. e.. 



[N, Vp] = [R, Vp] = [P, Vp] = [L, Vp] = 0. 



(144) 



Combining these results, we find the transformed Hamilto- 
nian in the accelerated and non-rotating frame as 

H"{t) = J d\"a\r")H'Jj^{r",p",t)a{r") + Vp. (145) 

All gauge contributions are contained in the definition of the 
modified single-particle Hamiltonian 

H'jy, p", t) = H'^p{M^v" + -r!, Af^p" + V\ t) 

+ S' +V'n' + V' (M^r") - n' (A/Tp") (146) 

+ P" (f^. xr"). 



By splitting the transformed Hamiltonian from Eq. J146> into 
constant, linear and higher order polynomial contributions in 
terms of the position f' and the momentum p', one finds 



ff^;(r", p", t) = [« (A/Tp") + f3 (A/^r") + 7] 

."2 ^ ^ (148) 

mVt fr",Af7e +SR,< 



P 

2m 



In there, we neglected terms of third and higher order in the 
gravitational potential. The coefficients read 



(149) 



13 =T' -v' xn 
d 



m-^[Vt{M{1l - p'),t) + VgiM{n +U),t)]. 

(150) 



and 



7-5' 



(151) 



If the constraints are satisfied identically, i. e., {a, (3, 7} — > 0, 
one obtains the single-particle Hamiltonian in the accelerated 
frame. 



JI2 



i7^;(r", p", t) + m Vt ( v\MTi + 5R, i ) , 



(152) 



provided the frame parameters as well as the center-of-mass 
of the drop capsule {"R. (t), "P {t), satisfy the classical 

equations of motion 



dL 



dp'' 



(153) 



sp 



V = = m(n' + ^i,x {n , V' = ^r- 
dn' ^ ^ dTZ' 

(154) 

With Eqs. ( I153t and jl54L we have recovered the Lagrangian 

L'[p', p' , iz' , TZ' , t] of Eq. (I122> . and one can express the ac- 
crued dynamical phase as an action integral 



2. Application of the frame transformation to the many-particle 
Hamiltonian in harmonic approximation 



to 



dt'L',p[n\n\t' 



to- 



(155) 



We will now apply the considerations of the previous sub- 
section to the falling trapped interacting many-particle system 
in the rotating frame. The single-particle Hamilton operator 
follows from the classical considerations of Eq. ( I124> and is 
given by 

H'sp{r',p',t) = ^ + Vt{M{t){r' - p'),t) 

+ Vg{M{t){f' + K), t) - p' {n, X f) . (147) 



V. TWO-COMPONENT ATOMIC GAS 

Let us now consider identical atoms with two internal elec- 
tronic degrees of freedom. In this section, we would like to 
demonstrate that the dynamics of these internal states decou- 
ples from the external motion of the trapped falling ultra-cold 
quantum gas, even in case of binary collisions, if the self- 
scattering and the cross-component scattering properties are 
identical. 
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A. Description in second quantization 

Henceforth, the inclusion of two degrees of freedom for the 
inner states in the definition of the field operators is required. 
These operators are denoted by a^{r), ji G {1, 2}, and satisfy 
the commutation relation 

[ap(ri), aj;,(r2)]^ = (5^,. 5{vi ~ V2). (156) 

The pair potential Vp, which is assumed to be local and ex- 
hibits translational and rotational invariance, can be written as 

Eh 



Iri-r^l) 



Mi..M4=i (157) 
«L (''i) ("^2) (i'2) ap4(ri)- 

The elements of \/a'ia'2A'3M4 linked to the self-species and 
cross-component scattering lengths of the atom in the dif- 
ferent states. In ^^Rb these quantities are - to a good ap- 
proximation - all equal, with a deviation of 3-4% fs?], so 
we can neglect the higher multipole contributions beyond the 

monopole (J=0) term V^;^^'-^'^^'^ « 5^,,^,,5^,,^,5^,^^Vy=''K 
This leaves us with writing Vp approximately as 

^P = \l d^E«]^('-i)«i('-2)^P°^(|ri-r2|)a.(r2)a^(ri) 

(158) 

We made use of Eq. ( I156> , and introduced the total particle 
density 

2 

n(r)^^at(r)a^(r). (159) 

We would like to point out that due to our assumptions Vp is 
now SU(2) invariant |53|. 

The two internal states can be coupled by a classical, 
traveling-wave laser field via electric dipole interaction \ 54\. 
The strength of the coupling is determined by the Rabi fre- 
quency which may be time-dependent if the laser is 
pulsed. A(t) denotes the, possibly time-dependent, detuning 
of the laser with respect to the resonant transition between the 
two levels. We want to consider large detunings in order to 
neglect any mechanical recoil effects of the laser on the atoms 
and, as usual, we have switched to an interaction picture os- 
cillating with the laser frequency [S^]. 

In order to describe the dipole interaction of identical parti- 
cles in the language of second quantization, we introduce the 
operators 



(160) 



with the well-known spin-l/2-matrices s^'^ i G {1, 2, 3}, Ap- 
pendixini Explicitly, the operators read 



S2 



dV 



S. = - / d^r 



4(i")ai(r) + al(r)a2(r) 
4(i")ai(r) - a\(Y)a^{Y) 
ai(r)ai(r) - a\{Y)a^(Y) 



(161) 
(162) 
(163) 



These operators fulfill the commutation relations of the angu- 
lar momentum operators for spin- 1/2 -particles and therefore 
represent the SU(2) symmetry. In other words, the proper- 
ties of the Pauli matrices translate to the picture of the second 
quantization. Clearly, the angular momentum algebra 



Si, Sj 



= iSijkSk (164) 
holds. All the operators R, P and L commute with Si, i.e. 

= 0. (165) 



S^,N 




Si, R 




S^,P 




Si, L 



Please note, that the definition of operators N, R, P, and L 
has to be modified, taking into account the two different inter- 
nal states, i.e. 



A^Y. [d'ral{Y)Aap{Y), 

. 1 n 



(166) 



/j=l,2 



where A e {N, R, P, L}, and A e {1, r, p, r x p}. 

We are now able to write down the full many-body Hamil- 
tonian, containing the one-particle, external dynamics (Hsp), 
the one-particle internal two-level dynamics (Vd), describing 
the interaction between matter and light, and the two-particle 
collisions (Vp), 



In there. 



H = [Hsp + Vd) + Vp. 



H^p^Y I d'ral{Y)H,papiY). 

, 1 n 



Hsp can be chosen as in Eq. J124t . Vd is given by 

Vd^Vdit) = ^it)Si + A{t)S3. 



(167) 



(168) 



(169) 



B. Separation of two-level dynamics and center-of-mass 
motion 

The dynamics of the many-particle system is determined by 
the Schrodinger equation 



(170) 



with H from Eq. ( I167> . As we have seen in the preceeding 
sections, a more efficient description of the dynamics can be 
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obtained by performing a frame transformation. In the present 
case, we have to account for both external and internal dynam- 
ics. Therefore, it would be favorable to separate both types of 
motion from each other We choose 

\^{t)) = U, Us^it) Up'it) Un'it) , (171) 

cf. Eq. ( fT39i . and 



canonical transformations in a time-dependent way will ulti- 
mately separate all "trivial" dynamics, including wave-packet 
spreading, from the essential many-particle physics. This is 
work in progress. We have also not touched questions of rel- 
ativity. This was done by intention in order to clarify all non- 
relativistic effects first, which by themselves are highly non- 
trivial and presumably dominant. 



Ud = T \ cxp 



dt'Vdit') 



(172) 



where T guarantees time ordering. The latter transformation 
formally eliminates the two-level dynamics, since it cancels 
the contribution Vd- The decoupling of the transformation Ud 
from the remaining frame transformation results from the fact 
that 



Vd, Hsp + Vp 







(173) 



due to Eq. ( I165t as well as the SU(2)-invariance of Vp, cf. 
preceeding subsection. In principle, the propagator Ud can be 
determined from the two linear ly indep endent solutions of the 
time-dependent Rabi problem I5j.l5a. 15711 . 

VI. CONCLUSIONS 

We have described the dynamics of an ultra-cold quantum 
gas in a long distance free-fall experiment. Starting from the 
classical mechanics of the drop capsule and a single trapped 
particle, we developed the quantum-field theoretical descrip- 
tion of a trapped, interacting degenerate quantum gas in a drop 
experiment in an inertial frame, the corotating frame of the 
Earth and the comoving frame of the drop capsule. By intro- 
ducing suitable coordinate transformations, it was possible to 
eliminate non-inertial forces and to focus on effects that take 
place on the mesoscopic length scale of the Bose or Fermi gas. 
The exact cancellation of non-inertial forces requires transla- 
tional invariance, the isotropy of the binary collisional poten- 
tial and the presence of a quadratic single-particle Hamilto- 
nian. This is well satisfied for ^^Rb, and the harmonic ap- 
proximation of the gravitational potential around the center- 
of mass of the BEC wave-packet is an excellent assumption. 
Corrections to it could be easily calculated perturbatively. 

If the atoms are two-level systems and coupled by an off- 
resonant traveling-wave laser field, this internal dynamics can 
be separated from the external motion, provided all scattering 
lengths are identical (SU(2)-invariance). 

This formalism provides us with an efficient way to de- 
scribe free-fall experiments, especially for numerical studies. 
It almost goes without saying that it is in particular valid and 
useful on the mean-field level. 

While we have discussed the Euclidean transformations 
corresponding to translation and rotation, we have omitted the 
scaling transformations |30], which are obviously very use- 
ful to model the adiabatic expansion of a BEC from a trap 
fsSi] , The consideration of the complete set of generalized 
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APPENDIX A: TAYLOR EXPANSION OF THE 
GRAVITATIONAL POTENTIAL 

The Taylor expansion of the gravitational potential 

V,{y' + p, t) . [1 + r' A + l-f + .. .] V,{p, t) 
= V,{p,t) + V^'\v',p,t) + V^^Hv',p,t) 



+5V^^\v',p,t). 

(Al) 

around the center-of-mass of the wavepacket defines a gradi- 
ent field, Eq. (I27> . a symmetric Hessian tensor, Eq. (I28> and 
the remainder forms a residual potential 



^V^f)(r',p) = ^FW(r',p), 



(A2) 



1=3 



which starts off with a leading third-order correction as 
Ff)(r',p) = l(r'®r'®r')-z;f (p), 



(A3) 



In order to obtain estimates for the Taylor coefficients at the 
surface of the Earth, we use the isotropic gravitational poten- 
tial and its derivatives 

Now, if this potential is expanded around a point p, we obtain 
the well-known multipole expansion in terms of scalar, dipolar 
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and quadrupolar components 

V,ir)=V,{W + p\)^V,ip)+r'^^ 



+ i(r'«.r') 



P dp 



,pdVg_ 



dp^ p dp 



(A5) 



,p® p 



1 



APPENDIX B: USEFUL RELATIONS FOR THE 
CONSTRUCTION OF THE LIE GROUPS ELEMENTS 

1. Euclidean space 

The angular momentum matrices (L — 1) ^^=1,2.3 satisfy 
the angular momentum algebra of Eq. ( I40l l. In a Cartesian ba- 
sis, the matrix representation is [ii\jk = —iT^ Qjk or explicitly 

£:L = h\ -i \ , e2^h\ 00|, (Bl) 
,0 i \-i 0) 



'0 -i 0^ 
i3 = h\ i I . 
,0 0; 



(B2) 



A finite orthogonal rotation is formed from these generators 
by 



r = r + Q X r 



(B3) 



For L = 1/2, the spin matrices s'^'-' satisfy the angular mo- 
mentum algebra of Eq. ( I40t as well. They are given by 



"2 " 2 U 



2 2 lO -1/ 



(B4) 
(B5) 
(B6) 



In there, tr'-*^ are the Pauli spin matrices 1 54.1 . 



2. Single-particle Hilbert space 

If the unitary operators act on the Heisenberg position and 
momentum operator, we find 

UlrUT^ = r + TZ, C/j, p C/^ = p + (B7) 

U^tUq^ Mq r, U^pUq^ Mq p. (B8) 

Please note that we have suppressed here extraneous hats to 
distinguish them from ordinary vectors in Euclidean space. 



3. Many-particle Fock space 

In the course of our calculations, we make use of some 
more basic commutator relations 



a{r),N =a(r), a(r),R =ra(r), (B9) 
a(r),pj = pa(r), [a(r), tj = 1 a(r). (BIO) 



They induce the unitary representations of the translational 
and rotational group 

^7^ ^ e^^/'^^f', = e-*/''2L, (b12) 

and we obtain again a representation of the group operations 

11^1111-^ = K + IZN, UIPU^^P + VN, (B13) 
C/^ R fjg = Mq R, uIPUq = Mq p. (B 14) 



APPENDIX C: CLASSICAL TRAJECTORY OF THE DROP 
CAPSULE IN A ROTATING FRAME 

In Section IV, the free-fall experiment in a rotating frame is 
discussed. A sketch of the physical situation is given in Fig.|5] 
6 and Q are linked via 



tan6' = — tanQ. 

a 



(CI) 



The classical trajectory of the drop capsule in the rotating 
frame of the Earth is given by 



Px (^) = 9 sin 6 cos ( 



t'^ 1 - 008(2175*) 

2" An? 



, , , sin 6 , sin(2il.t) , 




20. 



2r!. 



-' (,\ h 9.2^ - ■ 2a f^^ 1 - cos{2n,t) 
Pz{t) = h--t +5sm e[— 



(C2) 
(C3) 
(C4) 



where ft^ — and g = |g|. p'{t) solves Eq. il 15> . if the 
capsule is released with the initial conditions p'{t = 0) = 
(0, 0, h)^ , p'it = 0) 0, and if the centrifugal term is ne- 
glected. 

The rotation matrix Mg evolves around the axis fij with 
an angle Q — il^t (see Fig.|5}- Explicitly, it is given by 



cos Q cos 6* — sin Q cosQsin0> 
Mq — sin Q cos 6 cos Q sin Q sin 6 
— sin 61 cos 6 



(C5) 
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